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The atomic-scale spin structure of individual isolated skyrmions in an ultrathin film is investigated 
in real space by spin-polarized scanning tunneling microscopy. Their axial symmetry as well as their 
unique rotational sense is revealed by using both out-of-plane and in-plane sensitive tips. The size 
and shape of skyrmions change as a function of magnetic field. An analytical expression for the 
description of skyrmions is proposed and applied to connect the experimental data to the original 
theoretical model describing chiral skyrmions. Thereby, the relevant material parameters responsible 
for skyrmion formation can be obtained. 


Skyrmions are spatially localized solitonic magnetic 
whirls with axial symmetry and fixed rotational sense 
[1-4]. They have recently been observed in non- 
centrosymmetric bulk crystals [5-7] as well as in ultra- 
thin transition metal films on heavy-element substrates 
[8, 9], in which a sizable Dzyaloshinskii-Moriya interac¬ 
tion (DMI) [10, 11] induces their formation due to the 
lack of inversion symmetry. Magnetic skyrmions are in 
the focus of current research because they offer great 
potential as information carriers in future robust, high- 
density, and energy-efficient spintronic devices [12]: in 
addition to their protected topology and nano-scale size, 
they can easily be moved by lateral spin currents [12-16] 
and written as well as deleted by vertical spin-current 
injection [9]. In general, skyrmions couple very effi¬ 
ciently to spin currents and respond sensitively to spin 
transfer torques. For instance, skyrmion-based racetrack- 
type memory concepts [17] would profit from the fact 
that skyrmions can be moved with spin-current densi¬ 
ties being 5 to 6 orders of magnitude smaller than those 
needed to move magnetic domain walls [12]. The fabrica¬ 
tion of such skyrmionic devices using ultrathin films and 
multilayers can stay fully compatible with state-of-the- 
art technology. While numerous theoretical and simula¬ 
tion studies have concentrated on individual skyrmions 
and their physical properties, no high-resolution experi¬ 
mental characterization of the internal spin structure of 
skyrmions has been reported up to now, even though 
the knowledge about their actual size and shape provides 
the foundation for predictions about the interactions of 
skyrmions with spin currents or their manipulation by 
external fields as envisaged in potential skyrmion-based 
device concepts. 

In this letter we investigate the magnetic field depen¬ 
dent spin structure of isolated magnetic skyrmions by 
spin-polarized scanning tunneling microscopy (SP-STM) 
[18]. In contrast to techniques like magnetic force mi¬ 
croscopy [19] and electron holography [7], the spin struc¬ 
ture is measured directly rather than the magnetic field 
emerging from it. The combination of atomic-scale res¬ 
olution with spin sensitivity and the ability to mea¬ 
sure in external magnetic fields as high as several Tesla 
makes SP-STM an ideal tool to study the magnetic field- 


dependent spin structure of single nano-scale skyrmions. 
The evolution of their size and shape as a function of ex¬ 
ternal magnetic field is connected to micromagnetic the¬ 
ory of skyrmions [2, 3] via a proposed analytical descrip¬ 
tion of skyrmions. In this way, the values of the relevant 
material parameters can be extracted. 

Our experimental setup is sketched in Fig. 1(a): A 
magnetic probe tip with a well-defined spin orientation 
at the front atom is used in order to be sensitive to the 
spin-polarized contribution to the tunnel current, which 
depends on the projection of the local sample magneti¬ 
zation onto the quantization axis provided by the mag¬ 
netization direction of the tip. Here, we use an anti¬ 
ferromagnetic bulk Cr tip to avoid magnetic interactions 
of the tip with the sample or an applied field [21, 22]. 

As sample system we have chosen the bilayer of PdFe 
on an Ir(lll) single crystal substrate, which shows the 
typical magnetic field induced skyrmion lattice phase [9] . 
At the measurement temperature of T = 4.2 K the sam¬ 
ple exhibits pronounced hysteresis, enabling an investiga¬ 
tion of isolated skyrmions in a wide magnetic field range. 
Fig. 1(b) shows an overview of a PdFe area, exhibiting 
several circular skyrmions and two 360° domain wall sec¬ 
tions remaining from the spin spiral phase (top right). 
Due to the use of an out-of-plane sensitive SP-STM tip, 
the axisymmetric character of the skyrmions becomes di¬ 
rectly evident. 

When an in-plane sensitive SP-STM tip is used, the ap¬ 
pearance of the skyrmions changes: now two lobes with 
maximal and minimal spin-polarized current flow are im¬ 
aged per skyrmion (Fig. 1(c)) as a result of a positive or 
negative projection of the local magnetization direction 
of the skyrmion onto the spin direction of the tip. When 
tip and sample magnetization directions are orthogonal 
to each other, the spin-polarized contribution to the tun¬ 
nel current, and thus the magnetic signal, vanishes: for 
the image in Fig. 1(c) this is true close to the center of 
the Skyrmion and above and below the center. The two 
skyrmions in the sample area of Fig. 1(c) appear iden¬ 
tical, which is always the case for all skyrmions imaged 
with a given SP-STM tip. This implies that they ex¬ 
hibit indeed a unique rotational sense [23]. According to 
the symmetry selection rules of the DMI, these interface- 
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FIG. 1. Spin structure of individual skyrmion in 
PdFe/Ir(lll). (a) Sketch of the experimental setup of a spin- 
polarized STM tip probing a magnetic skyrmion. (b) Topo¬ 
graphic constant-current SP-STM image measured with out- 
of-plane sensitive magnetic tip; each blue circular entity is a 
skyrmion (U — +200mV, I — 1 nA, T — 2.2K, B — —1.5T). 
(c) Magnetic signal (methods [20]) of two skyrmions mea¬ 
sured with an in-plane magnetization of the tip, tut, re¬ 
vealing a two-lobe structure (U — +250 mV, I — InA, 
T — 4.2 K). (d) Same area as in (c) with inverted mag¬ 

netic field; due to the preserved rotational sense, the con¬ 
trast is inverted, (e), (f) Line profiles across a skyrmion 
along the rectangles in (c), (d), respectively, and fits with 
Eq. (1) ((e) c — (0.90 ± 0.01) nm, w = (1.18 ± 0.02) nm; (f) 
c = (0.91 ± 0.01) nm, w — (1.17 ± 0.01) nm) and correspond¬ 
ing calculated out-of-plane magnetization m z . The sketches 
show spins with atomic distance, colorized according to the 
SP-STM contrast. 


induced skyrmions are expected to be cycloidal (sketch 
in Fig. 1(a)) [11, 24], in agreement with recent density 
functional theory (DFT) calculations and Monte Carlo 
simulations for this system [25]. When the external mag¬ 
netic field which induces the skyrmions is applied in the 
opposite direction, the contrast of the two lobes of the 
skyrmions is inverted (Fig. 1(d)) since each spin in the 
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FIG. 2. Evolution of the size and shape of skyrmions 
in PdFe/Ir(lll) as a function of external magnetic field. 
(a)-(d) Magnetic signal of SP-STM differential conductance 
maps (methods [20]) with in-plane magnetized tip (U = 
+20 mV, I — 3nA, U mo d — +2.4 mV, T — 4.2 K) at the 
magnetic fields as indicated (Supplementary Movie 1 shows 
full data set [20]). (e) The size and shape of the skyrmion in¬ 
dicated by the box in (a) is evaluated by a fit with Eq. (2) as 
a function of magnetic field. Inset shows geometrical meaning 
of c, w and of d , which is numerically calculated; dashed blue 
line is a fit to d with 1 /(B — Bo). Solid black lines are ob¬ 
tained theoretically for the fitted set of material parameters 
A, D, and K (see text). Error bars correspond to standard 
deviation of fit parameters. 


sample is inverted while the spin structure of the anti¬ 
ferromagnetic tip remains unchanged, providing an ad¬ 
ditional proof for the unique rotational sense caused by 
the DMI. The two skyrmions in Fig. 1(c) and (d) appear 
at identical positions due to pinning at atomic defects. 

To characterize the size and the shape of a skyrmion 
we take height profiles across the center (see black rect¬ 
angles in Fig. l(c),(d)). Since there is no exact analytical 
expression to describe skyrmion profiles, we approximate 
the cross-section of a skyrmion using a standard 360° 
domain wall profile [26, 27], 
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where 0 defines the polar angle of the magnetization at 
position p, and c and re resemble the position and width 
of two overlapping 180° domain walls, respectively. To 
evaluate the measured data we include the projection of 
tip and sample magnetization in the fitting procedure 
(methods [20]). 

The agreement between experimental data and fit, see 
black data points and red fit lines in Fig. 1(e) ,(f), jus¬ 
tifies the chosen description. Furthermore, a compari¬ 
son to numerically calculated skyrmion profiles leads to 
the conclusion that Eq. (1) is an excellent approximation 
for a wide range of material parameters and field values 
(Fig. SI [20]). From Eq. (1) it is straightforward to deter¬ 
mine the perpendicular magnetization component m z (x), 
see blue dashed lines in Fig. l(e),(f), and the diameter of 
the skyrmion d, which we define as the diameter of the 
circle with m z = 0. Exploiting the axial symmetry, the 
spin structure of an isolated skyrmion in two dimensions 
is then described by: 


/-sin (6(p,c,w)) -x/p\ 

S{x,y)= -sin (6(p,c,w)) ■ y/p , (2) 

\ cos (9(p,c,w)) J 

where p = ^x 1 + y 2 is the radial distance from the center 
of the skyrmion located at the origin. Note that within 
this model the spin structure of the skyrmion is fully 
determined by only two parameters, c and w. 

The impact of external magnetic fields onto the size 
and shape of a skyrmion becomes evident in field- 
dependent SP-STM experiments: Fig. 2(a)-(d) show the 
identical sample area imaged at different external mag¬ 
netic field strengths. The color scale resembles the mag¬ 
netic contribution to maps of differential tunnel con¬ 
ductance (d//dC7), and the use of an in-plane sensi¬ 
tive tip leads again to the two-lobe appearance of the 
skyrmions. The decrease of the skyrmion size with in¬ 
creasing field can be directly seen in the displayed im¬ 
age sequence (Supplementary Movie 1 shows full data 
set [20]). For a quantitative analysis we fit a single iso¬ 
lated skyrmion (black box in Fig. 2(a)) with our two- 
dimensional skyrmion model and obtain the characteris¬ 
tic parameters c, re, and d. 

The evolution of these parameters with external field 
is shown in Fig. 2(e). The diameter of the skyrmion 
roughly scales with 1 /(B — Bq), see dashed blue line, 
in agreement with numerical calculations [3, 28]. While, 
in the investigated field range, the diameter changes by 
more than a factor of two, the effect on the width of 



FIG. 3. Validation of material parameters via micromagnetic 
simulations, (a) Internal spin structure of the skyrmion as de¬ 
scribed by Eq. (1) for the field values indicated. Inset shows 
visualization of spins with atomic distance as parametrized 
by Eq. (2); c, w given by fits in Fig. 2(e). (b) Comparison 
of experimental and simulated height profiles across an in¬ 
dividual skyrmion (box in insets) for several magnetic field 
values, and a fit with Eq. (1). Left and right insets show 
SP-STM experimental data from Fig. 2(a)-(d) and micromag¬ 
netic simulations based on the derived material parameters, 
respectively. 


the transition region is only about 25%. Consequently, 
as can be seen in Fig. 3(a), the skyrmion shape changes 
qualitatively with magnetic field, leading to a significant 
decrease in the number of spins with a component op¬ 
posite to the magnetic field. This results from a subtle 
balance of all involved energies, where the Zeeman en¬ 
ergy leads to a compressing force and the DMI stabilizes 
skyrmions against collapse to the ferromagnetic state. 

To assess these underlying interactions for the biatomic 
PdFe layer in the framework of micromagnetic continuum 
theory, we establish a connection to the standard energy 
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functional in cylindrical coordinates [1-3, 28-30]: 



We regard the precise knowledge of the spin structure of 
skyrmions as a prerequisite for further explorations of 
their possible application in spintronic devices. 

We thank A. N. Bogdanov, S. Rohart and J. Hage- 
meister for discussions. We acknowledge financial sup¬ 
port from the Deutsche Forschungsgemeinschaft (DFG) 
via grants SFB668 and GrK 1286. 


— K cos 2 0 — B z Ms cos 6 


pdp 


where exchange stiffness A, DMI constant D, uniaxial 
effective anisotropy constant K , and saturation magne¬ 
tization Ms are the material dependent parameters, B z 
is the external out-of-plane magnetic field and t is the 
film thickness. From DFT calculations [25], we estimate 
an Ms ~ l.lMAm -1 (methods [20]). The magnetiza¬ 
tion profile 6(p) is given by our experimentally verified 
skyrmion model, Eq. (1). Now for each set of A , D, iC, 
B the energy functional can be minimized with respect 
to c and w. These theoretical curves c(B) and w(B) are 
fitted to the experimentally obtained values for c(B) and 
w(B) via an error weighted least square fit with A, D, K 
as fitting parameters (Fig. 2(e)). The solid lines are the 
calculated values of d, c, w for A = (2.0 db 0.4) pj m -1 , 
D = (3.9 ± 0.2) mJ m -2 and K = (2.5 db 0.2) MJ m -3 
[31] as a function of magnetic field. These parameters 
are in the range expected for thin-film systems [25, 32], 
and the agreement of c and w obtained from theory with 
those from a fit to the experimental data is evident. 

To demonstrate that these derived material param¬ 
eters can be used to accurately reproduce the exper¬ 
imental data, we perform micromagnetic simulations 
[20, 30, 33]. Fig. 3(b) shows height profiles across an 
isolated skyrmion at four different magnetic field val¬ 
ues. The SP-STM data (black circles), the skyrmion fit 
with Eq. (1) (red line) and the height profile across the 
skyrmion in a micromagnetic simulation (blue dashed) 
nicely coincide, and the real-space agreement between 
experimental data and simulation is demonstrated in the 
insets to Fig. 3(b). Additionally, the accurate description 
of the field-dependent magnetism of the PdFe bilayer by 
the derived material parameters extends to lower mag¬ 
netic fields (Fig. S2 [20]). 

The presented SP-STM study provides access to the 
actual spin structure of an isolated skyrmion, enabling 
a direct comparison to micromagnetic theory. Since 
the field-dependent evolution of size and shape of single 
skyrmions is governed by a balance of the magnetic in¬ 
teractions, a detailed experimental characterization can 
yield the relevant material parameters such as the DMI, 
which is responsible for the stability of these particle¬ 
like states. This procedure is applicable to all DMI- 
stabilized magnetic objects; however, for thick layers or 
small DMI the demagnetizing field should be taken into 
account explicitly to obtain realistic material parameters. 
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